Nonlinear Evolution of Anisotropic Cosmological Power 
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There has been growing interest in the possibility ol testing more precisely the assumption of 
statistical isotropy of primordial density perturbations. If it is to be tested with galaxy surveys at 
distance scales < 10 Mpc, then nonlinear evolution of anisotropic power must be understood. To 
this end, we calculate the angular dependence of the power spectrum to third order in perturbation 
theory for a primordial power spectrum with a quadrupole dependence on the wavevector direction. 
Our results suggest that primordial power anisotropies will be suppressed by < 7% in the quasilinear 
regime. We also show that the skewness in the statistically anisotropic theory differs by no more 
than 1% from that in the isotropic theory. 

PACS numbers: 98.80.-k 



It is commonly assumed that primordial density per- 
turbations are statistically isotropic, and statistical 
isotropy is a prediction of most structure-formation theo- 
ries. The notion of statistical isotropy can be quantified, 
though, by considering models in which it is broken, and 
a growing literature has discussed physical models that 
produce primordial perturbations that are not statisti- 
cally isotropic [l[. 

The manifestations of departures from statistical 
isotropy can take on many forms, but one simple ex- 
ample Q predicts a primordial power spectrum P(k) 
with a quadrupole dependence on the direction k of the 
wavevector fc, 



P(k) = A(k) [1 + gM/ik)], 



(1) 



where V-z{x) = (3x 2 — l)/2 is the second Legendre poly- 
nomial, fik — k ■ e, and e is a preferred direction. 1 Refer- 
ence Q constructed a CMB minimum- variance estimator 
for the power-anisotropy amplitude g* and showed that 
the Planck satellite will be sensitive to a value of <?* as 
small as ^0.02, a number which can be estimated roughly 

— 1/2 

by ~10iV pix ' , where 7V pix is the number of statistically 
independent pixels on the sky detected by Planck. 

Similar arguments suggest that the sensitivity of a 
galaxy survey, like the Sloan Digital Sky Survey [J] or 
Two-Degree Field [H , should have a comparable sensitiv- 
ity. One issue that will arise, though, in testing statisti- 
cal isotropy of primordial perturbations is the quasilinear 
evolution of the power spectrum. The root-mean-square 
density-perturbation amplitude becomes of order unity 
at distance scales ~10/i _1 Mpc, and so quasilinear evo- 
lution of the density field must be taken into account 
if the mass distribution measured on these scales in the 
Universe today is to be used to infer the primordial mass 
distribution. 



In this paper, we study the nonlinear evolution of den- 
sity perturbations to see how nonlinearity affects statisti- 
cal isotropy. Does nonlinear evolution amplify departures 
from statistical isotropy? Or possibly suppress them? To 
take the first steps to address this question, we calculate 
the power spectrum, to third order in perturbation the- 
ory, under the assumption that primordial perturbations 
have the form given in Eq. ([I]). The bottom line is that 
quasilinear evolution suppresses departures from statisti- 
cal isotropy, but by only < 7% compared with the linear 
theory. Thus, galaxy surveys in the linear/quasilinear 
regime will still be useful diagnostics for departures from 
statistical isotropy. We also calculate skewness, finding 
that the quadrupole power anisotropy changes it by no 
more than 1%. 

To proceed, we use third-order perturbation theory 
to determine whether a primordial power anisotropy is 
amplified or suppressed in the quasilinear regime. The 
power spectrum P(k, z) is defined by the ensemble aver- 
age of the two-point correlation of the Fourier transform 
<5(fc, z) of the density perturbation through 



(S(k 1 ,z)S(k 2 ,z)) 



(2TT) 3 S D {k 1 



fca)P(fc x 



(2) 



where Sp is the Dirac delta function. The density per- 
turbations can be expanded in terms of the linear-theory 
density-perturbation amplitude, which has a redshift de- 
pendence proportional to the linear growth factor D(z), 



6(k,z) = J25 n (k)D n (2 



(3) 



Note that our g* is 3/2 times that in Ref. 0. 



In the linear regime, each Fourier mode grows at the same 
rate, and so the linear-theory power spectrum Pn n (k), 
corresponding to the linear density perturbation #i(fe), 
evolves in such a way that the quadrupole dependence 
of the primordial power spectrum [Eq. {!])] is preserved. 
However, the nonlinear power spectrum will have differ- 
ent anisotropy structure, and we shall investigate it fol- 
lowing prescription developed in Refs. (y, LZI, [8|. 

Reference |6| showed that solution of the nonrelativis- 
tic fluid equations (i.e., continuity, Euler, and Poisson 
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equations) allows one to write the nth-order density- 
perturbation amplitude S n in terms of the linear pertur- 
bation Si through 



(2tt) 3 J (2tt) 3 
x {2TT) 3 F n {q 1 ,---,q n )5 1 (q 1 )- 



■5i(q n ). (4) 



For any n, the function F n can be obtained by recursive 
relations given in Ref. [6|. In particular, the expressions 
for F 2 and P3 that are directly relevant for our purpose 
are explicitly given by Eqs. (A2) and (A3) of Ref. [6J. 2 
As all odd moments of 5\(k) vanish, the power spectrum 
to third order is given as 

(2TT) 3 S D (k 1 +k 2 )P^(k 1 ,z) 

= D 2 (z)(S 1 (k 1 )S 1 (k 2 )) + D\z) [(5 2 (k 1 )5 2 (k 2 )) 
+<5 1 (fci)fc(fca)> + (<5 3 (fci)5i(fc 2 ))] . (5) 

The second term, which evolves as D 4 (z), is the next- 
to-leading term evaluated to the third order of density 
perturbation £3, while the first term is the linear part. 
We further define a quantity P mn as 

(5rn(kl)S n (k 2 )) = {2^f5 D {k 1 + fe 2 )P mn (fel), (6) 

and with this definition, Eq. ([5]) is rewritten as 

P^(k,z) ^ D 2 (z)P 11 (k)+D 4 (z)[P 22 (k)+2P 13 (k)]. (7) 

We then need to write P 22 and Pi 3 in terms of the 
linear spectrum P\ in — P\\. This procedure is straight- 
forward,_and the results are very similar to those given 
in Ref. 



P 22 (fc) 



dq 

(2^f 



F^(q,k 



x Piin(q)Piin(A: - q), (8) 
Pis(fe) = 3P lin (fc) J ^P 3 (s) ( < 7,- 9 ,fc)P lin (q), (9) 

where the symmetrized function P„ is obtained by sum- 
ming n\ permutations of F n over its n arguments and di- 
viding by n\. The only difference from the expressions in 
Ref. [8[ is that the linear power spectrum now depends 
on both magnitude k and direction k of wavevector k. 

In order to illuminate the anisotropy structure, we ex- 
pand Eqs. |8]| and (J9j) with Legendre polynomials. The 
most general power spectrum today (D = 1) can be ex- 
panded in Legendre polynomials V n (x) as 

00 

9* Pn 

(Vk)B mn (k), (10) 

m,n— 



2 Strictly speaking, F2 in Ref. [f| is for an Einstein-de Sitter Uni- 
verse, but Ref. [S|] shows that its form in a ACDM model with 
Q m ~ 0.3 differs by less than 1%; the same should be true for 
*3. 



where the sum on n is restricted to even positive inte- 
gers. To the next-to-leading order, the expansion can be 
written, 

P^(k) = A Un (k) + B oa (k) + g 2 B 20 (k) 

+ {gMUk) + B 12 {k)} + g 2 B 22 {k)} P 2 (fc) 
+ glB 2i (k)V i (k) 1 (11) 

where An n , defined by Eq. (frj, is the isotropic linear 
power spectrum. Expressions for the expansion coeffi- 
cients B mn (k) are given at the end of the paper. We 
then rewrite Eq. (fTTjl as 



P^(fc) = A^(k) [l + 9 f (k)V 2 (Hk) 

9 { 2 \k) = g*ci(k) + g 2 c 2 (k), 
gi 2) (k) = glc 3 (k), 



(12) 

(13) 
(14) 



with the following definitions of A^ 2 \ ci, c 2 , and C3 
A^(k) 

A^{k)a(k) 



A ria (k) + B 00 (k) + g 2 B 2Q (k), (15) 

" A Un (k) + B 12 (k) [i = l] 

B 22 (k) [» = 2] . (16) 
B 24 (k) [i = 3] 



We note that the primordial quadrupolc anisotropy also 
affects the isotropic part of the next-to-leading order 
power spectrum, A^ 2 \k). While this is interesting, that 
correction is expected to be very small, being suppressed 
by g 2 - Thus, hereafter we neglect this term; keeping it 



-.(2) 



and 



only gives correction of the order of g 3 and g 4 to g 2 
g±, respectively. 

For the linear power spectrum An n (k), we use the fit- 
ting formula for the transfer function given in Ref. [Io| 
with current values for the relevant cosmological param- 



eters 

If 5* <C 1, we may neglect all the terms propor- 
tional to g 2 , and thus c\(k) is the only quantity relevant 
for anisotropy structure. As linear theory simply gives 
c\{k) = 1, the quantity c\{k) — 1 represents enhancement 
of the quadrupole anisotropy due to quasilinear evolu- 
tion. In Fig.[U we plot c\{k) — 1. This figure shows that 
at large scales k < 10 -2 Mpc^ 1 , the nonlinear effect is 
subdominant and thus anisotropic structure is the same 
as that for the linear theory: c\(k) « 1. The quadrupole 
anisotropy then decreases at smaller scales and it reaches 
minimum (ci — 1 ~ —0.07) in the quasilinear regime, 
k ~ 0.1 Mpc" 1 . 

Third-order perturbation theory becomes less accurate 
for even higher wavenumbers, k > 0.1 Mpc -1 , at z = 0. 
However, it remains very accurate at such (comoving) 
scales in the high-redshift Universe, z > 1 [12J. For gen- 
eral redshifts z, the anisotropic power spectrum is given 
by Eqs. (fT!?)) - (flT)|) with replacements An n — > D 2 i K z)A\\ n 
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k [Mpc- 1 ] 



FIG. 1: Coefficients of anisotropic terms in the third-order 
power spectrum, ci(fc) — 1, C2(fc), and c$(k). The definitions 
are given in Eqs. (I12|| - (|14p . Solid curves are the result of nu- 
merical integration, while dotted curves are the fitting func- 
tions (fTTT) - (fT9|) . 



and B rnn — > D i (z)B mn . The enhancement of anisotropy 
is then given by c\(k,z) — 1, which is well approxi- 
mated by [ci(k) — l]D 2 (z) in the quasilinear regime where 
A\i n > \B mn \. It is -0.03 at k = 0.1 Mpc" 1 and z = 1; 
the suppression is weaker. One can also see amplified 
anisotropy for k > 0.3 Mpc -1 at redshifts when the third- 
order approach is still valid at such scales. 

The corrections on the order of g 2 are represented by 
C2(/c) and c^{k), both of which are also shown in Fig. [1] 
Besides it is suppressed by additional C2(k) is intrin- 
sically smaller than c\ — 1, thus giving only a minor cor- 
rection. While possessing different anisotropy structure, 
the higher multipole term proportional to c 3 {k) would 
also be small in the quasilinear regime. 

The results of our computations can be approximated 

by 

ci(fc) = 1 - 0.463A 2 (fc) + 0.886A 4 (fc) - 0.407A 6 (^,7) 
c 2 (fc) = -0.133A 2 (fc) +0.195A 4 (fc) - 0.079A 6 (fc),(18) 
c 3 (k) = -0.163A 2 (fc) + 1.062A 4 (fc) - 1.082A 6 (fc) 

+ 0.373A 8 (fc), (19) 

in terms of the linear-theory density-perturbation ampli- 
tude at wavenumber k, A 2 (k) = (k 3 /2ir 2 )Au n (k). These 
fitting functions are plotted as dotted curves in Fig. Q] 
They were obtained explicitly for the current best-fit cos- 
mological parameters 11], but when written in terms of 
A 2 (fc), should also be accurate for other cosmological pa- 



rameters. Strictly speaking, the validity of third-order 
perturbation theory breaks down when the terms higher 
order in A 2 (fc) become important; in this case, one will 
need to go to fourth-order or even higher-order correc- 
tions to obtain the evolved power spectrum. 

We now compute the skewness S3 = {^(5 3 ) / (<5 2 ) from 
the anisotropic primordial power spectrum, to second or- 
der in perturbation theory. The standard result, for an 
Einstein-de Sitter Universe, is S 3 = 34/7 [HEi], 

and the 

result for a ACDM Universe differs by less than 1%. The 
three-point correlation function at zero lag is a Fourier 
transform of (<5(fei, z)6(k2, z)8(k 3 , z)), and thus the lead- 
ing contribution comes from Z? 4 (z)((5i(fci)<5i(fc2)<52(fc3)), 
etc. Rewriting 82 in terms of the linear fluctuation #1 
and the integration kernel F% , and using the definition of 
linear power spectrum, we obtain, 



<<5 3 (z)) = 6D 4 (z) 



dki 



(2tt) 3 J (2tt) 

X fW(fcl)Ain(fc2 



J (fcl,fc 2 ) 



(20) 



After integrating over the directions of wavevectors, k\ 
and ^2, we find 



/' (y-yi 



dkk 2 

~2~^ 2 ~ 



Al in (fc) 



(21) 



On the other hand, the variance (5 2 ) can be computed 
in a similar manner, and to the same order, we obtain 

{5 2 {z))=D 2 {z) J Jt_ Plin (k)=D 2 (z) J ^-AUk)- 

(22) 

Therefore, we find the skewness to be 



5*3 



(6 2 (z)) 2 



34 

y 



175' 



(23) 



The requirement that the power spectrum be positive 
definite restricts the value < g* < 1, and so the largest 
possible deviation from the isotropic value of 34/7 is less 
than 1%. Thus, the skewness will not be an effective 
discriminant for anisotropic power. 

To conclude, we used third-order perturbation the- 
ory to calculate the quasilinear evolution of a primordial 
power spectrum. Our results show that nonlinear evo- 
lution suppresses the primordial anisotropy. The sup- 
pression is ^7% at k = 0.1 Mpc -1 and z = compared 
with the linear theory. This must be taken into account 
when interpreting the result of searches in the quasilin- 
ear regime for primordial anisotropy, and it indicates that 
nonlinear scales arc still valuable as tests for the isotropy 
of primordial power, since the suppression is weak. We 
also found that the standard prediction for the skewness 
is little changed if the primordial power spectrum has a 
quadrupole anisotropy. 

Our perturbative results seem to indicate that a 
quadrupole anisotropy is enhanced at smaller scales, but 
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this enhancement cannot be trusted at low redshifts, as it 
occurs at scales where our perturbative approach breaks 
down. It seems counterintuitive to think that growth of 
perturbations in the highly nonlinear regime can act to 
enhance the primordial anisotropy, but it will require an 
N-body simulation to be sure. Likewise, we imagine that 
departures from statistical isotropy that are higher order 
in angle (octupolc, etc.) will also be suppressed by quasi- 
linear evolution, but we leave that calculation for future 
work. 

Before closing, we provide explicit expressions for the 
quantities in Eq. (jTUJ) . For the numerical evaluation of 
Eqs. ijHJ and ([9|), we shall choose fc along the z-direction, 
fc = (0,0,1), and the preferred direction e in the x-z 
plane, e = [y/l — m| , 0, /ifc). We then use the spheri- 
cal coordinate for q = q{-\fl — /i 2 cos <fi, \J\ — /i 2 sin 4>, f£), 
and dq = q 2 dqdnd<f>. With this notation, B mn in Eq. (fTO]) 
are given as 



Sm.n. — 



inn 
A 



271* 



dq / djj. 



[7kn + q(3 - lOfi 2 )} 2 



(fc 2 



X A lin (q)Ay ln Ufk 

k 2 A Un (k) 



q 2 — 2kqfi) 3 
1 - 2kq») , 



dqA lin (q) 
k 



<J •r 



„2\3 



hi 



(26) 



\k-q 

where S mn and T mn correspond to P22 and P13, respec 
tively. Explicit forms of T mn , J mn , and IC mn are 



-*-00 

I12 
I22 

2^24 

J12 
J20 



fc 2 + q 2 - 2kq^ 
392 




(27) 


2q 2 (3/i 2 - 1) + (fc 2 - 


2kqn){3n 2 + 1) 


(28) 


784 




2q 2 -Akqfi + k 2 (3fj, 2 


-1) 


(29) 


3920 


J 


(fc 2 + q 2 )(3fi 2 - 1) - 


kqfiiSfi 2 + 1) 


(30) 


2744 
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54880 
— 8fc<7/i(5/i 2 

1 



[4fc 2 (3/? - 1) + g 2 (35/ - 30^ 2 + 3) 



1008 
1 

40320 



q 2 
fc 4 

9 V 



-3)] 
79 4 



2 4 



fc 2 

375 ^ 



5006+ 1732-^ 



9 4 
fc 4 



q 6 
fc 6 



(31) 
(32) 

(33) 



10 



•J22 



7 

18 



J24 



1 



fc 4 

, 90— + 

201600 V q 4 

q 4 q 6 

-540fr - 315f« 
fc 4 fc 6 



fc 2 
135— 

q 2 



1846 



268 h 



IC12 
IC20 



7q 2 + 2fc 2 
672fc 5 g 3 ' 

6fc 6 + 41fcV + 76fc 2 g 4 + 21g 6 

5376fcV 
7 7 
10 18 

6fc 6 + 25fc 4 g 2 + 20fc 2 g 4 + 21g 6 
26880fcV 



(35) 
(36) 

(37) 



and all the other components vanish. Note that T 00 is 
the same as Eq. (19) of Ref. 
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